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Abstract 

We derive a small-time expansion for out-of-the-money call options under an exponential Levy model, using 
the small-time expansion for the distribution function given in Figueroa-LopezfcHoudre |FLH09] . combined with a 
change of numeraire via the Esscher transform. In particular, we find that the effect of a non-zero volatility a of 
the Gaussian component of the driving Levy process is to increase the call price by \aH'^e^v{k){l + o(l)) as f -)■ 0, 
where v is the Levy density. Using the small-time expansion for call options, we then derive a small-time expansion 

for the implied volatility ot (k) at log-moneyness k, which sharpens the first order estimate at (k) ~ tiog(i/t) gi'^sn 
in |TnkvlO| . Our numerical results show that the second order approximation can significantly outperform the first 
order approximation. Our results are also extended to a class of time-changed Levy models. We also consider a 
small-time, small log-moneyness regime for the CGMY model, and apply this approach to the small-time pricing of 
at-the-money call options; we show that for Y G (1, 2), limt_>o t^^''^ E(St — 5*0)+ = S'oE*(Z+) and the corresponding 
at-the-money implied volatility at{0) satisfies limt-to at{0) /t^^^ — \/27rE*(Z+), where Z is & symmetric Y- 
stable random variable under P* and Y is the usual parameter for the CGMY model appearing in the Levy density 
v{x) = C2;-'-'^e-*^"l{,>o} +C|2;ri-'^e-'^i"ll{,<0} of the process. 

1 Introduction 

Levy processes have played an important role in the development of financial models which can accurately approximate 
the so-called stylized features of historical asset prices and option prices. In the "statistical world" , financial asset prices 
exhibit distributions with heavy tails and high kurtosis as well as other dynamical features such as volatility clustering 
and leverage. In the "risk-neutral world" , market prices of vanilla options exhibit "skewed" implied volatilities (relative 
to changes in the strike), contradicting the classical Black-Scholes model which predicts a flat implied volatility smile. 
The smile phenomenon has been more pronounced after the 1987 market crash. Concretely, out-of-the-money equity 
put options typically bear a higher risk-premium (larger implied volatilities) than in-the-money puts. This effect is 
more dramatic as the time-to-maturity decreases. As explained in |CT04j (see Section 1.2.2), the latter empirical fact 
is viewed by many as a clear indication that a jump risk is recognized by the participants in the option market, and 
stochastic volatility models are, in general, not able to reproduce the pronounced implied volatility skew of short-term 
option prices unless the "volatility of volatility" is forced to take high values. 

The literature on small-time asymptotics for option prices and implied volatilities has grown signiflcantly during 
the last decade. For recent accounts of the subject in the case of stochastic volatility models, we refer the reader 
to |GIILOW09| for local volatility models, jFJLlO] for the Heston model, |Forde09) for a general uncorrelated local- 
stochastic volatility model and [FordelO] for SABR type models. We concentrate here on asset price models with 
jumps. For an Ito semimartingale model for the underlying price process {St), CarrfcWu |CW03| argued, by partially 
heuristic arguments, that the price of an out-of-the-money call option converges to zero at sharply different speeds 
depending on whether the underlying asset price process is purely continuous, purely discontinuous, or a combination 
of both. For instance, in the presence of jumps, they argue that the 

nSt-K)+ -{S^-K)+^ct, (1) 
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for some constant c 7^ 0, as the time-to-maturity t tends to oQ while the call price converges at the rate 0(6""^/*) for a 
purely-continuous model. These statements were subsequently exploited in |CW03j to investigate which kind of model 
is more adequate to describe the observed market option prices near to expiration. They concluded the necessity of 
both a continuous and a jump component to describe the implied volatility of S&P 500 index options and argued, 
based on simulation experiments, that the theoretical asymptotic behavior is usually manifested by options maturing 
within 20 days. We also refer the reader to jAS02| for further empirical evidence on the presence of both a continuous 
and jump component. 

Using the closed- form expressions for call option prices, BovarchenkofcLevendorksii }BL02| (see also |Lev04a) . 
|Lev04b| ■ jLev04c| ) establish the following small-time asymptotic behaviour 

^¥.{St- K)+ ^ t I [e'' ~e'')+v{dx), (fc > & t ^ 0) , (2) 



5o 

for several popular exponential Levy models St = Sqc^* , where k is the log-moneyness k := log(7^/ Sq) and v is the Levy 
measure of the underlying Levy process (Xt). Subsequently, Levendorskii |Lev08| obtained ([2]) under certain technical 
conditions (see Theorem 2.1 therein), namely that / {\x\'^ Al)e^'^^^^^v{dx) < 00 for some e > 0, and limf_i.o E(<5't — 
exists in the "out-of-the-money region". More recently. Roper |RoplO| and Tankov |TnkvlO| prove that ([2]) holds for a 
general Levy process (Xt) under mild conditions, using the first-order small-time moment asymptotic result 

lim^E{^{Xt)}^ J ^{xHdx), (3) 

valid for functions (p that converges to as a; at an appropriate rate (see, e.g. Figueroa-Lopez jFL08| for details). 
In particular, it suffices that e^v{dx) < 00. |Lev08] also provides a natural generalization of ^ for a wide class 

of multi-factor Levy and Markov models. 

As a corollary of ©, |Ropl0| and |TnkvlO| prove independently that the implied volatility at{k) for exponential 
Levy models explodes near expiration for out-of-the-money vanilla options. This is a very peculiar feature of financial 
models with jumps (sec Remark 12.61 below for a brief discussion about its meaning). jTnkvlO] goes one step further 
and shows that 

as t 0. For at-the-money call option prices, |RoplO| also shows that the leading order term is 0{\/i) and docs not 
depend on the jump component of the model. Moreover, the at-the-money implied volatility converges to the volatility 
of the Gaussian component of the driving Levy process, and the limit is zero if the Levy process has no Gaussian part. 
For bounded variation Levy processes and for certain tempered-stable like Levy processes, [TnkvlO] gives also the 
first-order asymptotic behavior of at-thc-moncy implied volatilities. The asymptotic behavior (j4]) is in sharp contrast 
with a pure-continuous stochastic volatility model, where the implied volatility converges to a non-negative constant 
which depends on the shortest distance from zero to the vertical line with x-coordinate equal to the log-moneyness of 
the call option, under the Riemmanian metric induced by the diffusion coefficient for the model (see, e.g. |GHLOW09) , 
|FJL10| . }Forde09) . [FordelOj l. 

In this article, we extend previous results by computing the second order correction term ai(fc) in the call option 
price approximation: 

l^nSt-K)+ = ao{k) + ai{k)t + oit) {t ^ 0) . (5) 

t bo 

An important component in our proofs is played by the recent higher order small-time expansions for the distribution 
function of a Levy process obtained in Figucroa-L6pczfcIIoudre [FLH09j . In the spirit of the Black-Scholcs formula 
and the classical change of numeraire, our approach exploits an appealing representation of the prices of out-of-the 
money options in terms of the tail distribution functions of the underlying Levy process under both the original risk- 
neutral probability measure P and under the martingale probability measure P* obtained when we take the stock as 
the numeraire; i.e. F*{A) := E(S'(1^) (see e.g. Chapter 26 in |Bjr09| and references therein). The latter measure P* 
is sometimes called Share measure (see e.g CarrfcMadan jCM09| ) . Our results allow us to quantify precisely the effects 



^Actually, ICW03| wrote Ilt{K) — {So — K)+ = 0(t), even though in their empirical analysis they are assuming a stronger statement 
such as |(TJ. 
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of a non-zero Gaussian-component in the call option prices near expiration. Wc find that a continuous-component 
volatility of a will result in an call price increase of ^a^t^ef^viW) (per each dollar of the underlying spot price), where 
V is the Levy density and k is the log-moneyness. 

We also derive the corresponding small-time asymptotic behaviour for implied volatility, showing precisely how the 
implied volatility diverges to oo (see Section [2]) . We find that the dimensionlcss implied variance does tend to zero 
as we would expect, but very slowly; in fact slower than for any p > 0, and consequently the implied volatility 
explodes in the small-time limit. Furthermore, we characterize the asymptotic behavior of the relative error of the first 
order approximation, which is then used to obtain a second-order approximation for the implied volatility of out-of-the 
money call options. According to our numerical results (see Section [B] for the details), the second order approximation 
significantly reduces the error compared to that of the first order approximation, achieving up to a two-fold relative 
error reduction in some cases. 

We later extend our analysis to the case of a time-changed exponential Levy model Zt — with an independent 
absolutely-continuous time change Tj = J^Ysds satisfying some mild moment conditions (see Section [3]). The time- 
changed Levy model was proposed in |CGMY0"3] to incorporate the volatility clustering and leverage effects commonly 
exhibited by financial price processes. We show that the small-time behavior of call option prices depends not only on 
the triplet of the underlying Levy process X but also on the time-zero first and second moments of the speed process 
(Yt) and the quantity 

which is assumed to exist. In some sense, 7 measures the current average acceleration of the random clock. Under 
mild conditions, we show that 

\^¥.{St-K)+ = Eyoao(fc) + a^{k) + -fao{k)\t + o{t), (t 0) , 

where ao,ai are the first and second order terms appearing in the pure-Levy option price approximation ([5]). For a 
Cox-IngersoU-Ross (CIR) speed process 

dYt = n{e-Yt)dt + a^tdWu Yo=yo 

the current acceleration of the process is 7 = k{0 — y^) and, hence, call option prices will exhibit the following 
small-maturity asymptotic behavior: 

\ ^E(5t - K)+ = yoao(fc) + [ylai{k) + ^{9 - yo)ao(fc)] t + o(t), (t ^ 0) . 

As seen from this expression, a mean reversion speed of k will increase (resp., decrease) the call option price when the 
current volatility j/o is above (resp. below) the long-run mean volatility value 9. 

In Section m we also consider a small-time, small log-moneyness regime for the CGMY model of jCGMY0"2| . The 
CGMY model is a particular case of the more general KoBoL class of models, named after the authors Kop95| (who 
first introduced the symmetric version of the model under the name of "truncated Levy fiights" ) and |BL02| . Using 
the fact that {Xtlt^l^)^ converges weakly to a symmetric alpha-stable distribution with a = y as t — ^ 0, we show that 

limri/^^E(5t-^o)+=E*(^+), 

for y G (1,2), where Z is a symmetric y-stable random variable under P*. Wc then apply this result to small-time 
pricing of at-the-money call options for the CGMY model. Our method of proof is new and based on the following 
representation by Garr&:Madan |GM09j 

^nSt-K)+^r{Xt-E>\og^), (6) 

where E is an independent exponential random variable under P* with parameter 1. As a corollary, we conclude 
that the corresponding at-the-money implied volatility (Tt(0) satisfies lim(_i.o CTt(0)/t^/^^ 2 — v^27rE*(Z+). |TnkvlO| 
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obtains a similar result in a more general model using a different approach based on a Fourier-type representation for 
call option prices. Let us also remark that the method of proof introduced here can be applied to a large class of Levy 
processes whose Levy densities are symmetric and dominated by stable Levy densities, and behave like a symmetric 
F-stable process in the small-time limit (see Remark 14.31 for the details). 

In section [Sj we derive a similar small-time estimate for variance call options using the well known fact that the 
quadratic variation [X]t of a Levy process is itself a Levy process. Using the main resuh in [FL08] . we find that an 
out-of-the- money variance call option which pays {[X]t — K)^ at time t is worth the same as a European-style contract 
paying (In ^)'^ — K)+ at time t as t — > 0, irrespective of the Levy measure i^(-). The diffusion component of (Xt) does 
not show up at leading order for small t. See also jKRMKlT] for a related discussion on the difference between the 
small-time behaviour of variance call options on the exact quadratic variation and its discretely sampled approximation 
for Levy driven models. 



2 Small-time asymptotics for exponential Levy models 

Consider an exponential Levy model for a stock price process 

St = 5*06^* (7) 

where (Xt) is a Levy process defined on a complete probability space (f2,P, J^) with generating triplet ((t^,&, r^). We 
are assuming zero interest rate and dividend yield for simplicity!^ and that P represents a risk- neutral pricing measure. 
We assume that Jj^i^^ e''v{dx) < oo and that the following condition is satisfied 

^+9^'+/ (e" - 1 - 2:l|,|<i)i.(dx) = 0, (8) 

^ J -co 

so that St = Sqc^* is indeed a P-martingale relative to its own filtration. 

Throughout the paper, we also assume that the Levy measure i'[dx) admits a positive density, denoted J^(x), and that 
this density is in M\{0} satisfying supj^-i^^ v{x) < oo for every e > 0. The choice between the Levy measure ^{dx) 
and density I'ix) should be clear from the context. Under the previous standing condition, Figueroa-Lopez&Houdre 
|FLH09| show the following result (see Remark 3.3 and Proposition 3.4 therein): 

Theorem 2.1 \FLH09f . Let y > be fixed. Then, we have the following small-time behaviour for the distribution 
function of Xt : 

i >y)^ iy[y, c^) + ^td2{y) + ait) {t ^ 0) , 

where 

d2{y) = d2{y]b,<7,v) = ~a'^u'{y)+2bv{y)-v[y,'X)f + u{^y,yf 

-kv fv r rhy rv 

/ v{u)v{x)dudx — 2v{y) I xv{x)dx + 2v{y) / / [v{u) — v{y))i>(x)dudx . (9) 

oo J y-x Ji'y<|a;|<l J-^yJy-x 

Furthermore, if the pure-jump component of (Xt) has finite variation, i.e. J^^^^i \x\i'{dx) < oo, then d2 simplifies to 

.2,//„,\ , of,_,,/„A ,,r„, ,„\2 



d2{y) = d2{y',b,a,v) ^ —a v'{y) +2b{)v{y) — v[y,co) + / v{u)v{x)dudx — 2 I I i'{u)v{x)dudx , (10) 

Jy—x Jy J —oo 

where bo is the drift of the pure-jump component of (Xt) defined by bo = b — J^j.^<:i xv(dx). 



^For a non-zero constant interest rate r and dividend rate g, the results in this paper will not be qualitatively any different, because we 
can just replace the stock price process (St) with the forward price process (e~(''~'''*5t)t, which is a martingale (see, e.g. Chapter 11 in 
|CT04n . 
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Remark 2.1 The double integrals in ([9]) and (fTO|) are well-defined. For instance, by the symmetry of s(u)s{x) about 
the line u ~ x, 

V ry ry/'^ ry ry ry/'^ ry ry 

I v{u)v{x)dudx ~ / v{u)v{x)dudx + / v{u)v{x)dudx + / v{u)v{x)dudx 

Jy~x Jo Jy-x Jy/2Jy-x Jy/2Jy/2 

/•y/2 ry ry ry 

— "2 / i'{u)v{x)dudx + / v{u)i>{x)dudx (11) 

Jo Jy-x Jy/2Jy/2 

ry/'^ 

<2[ sup v{u)\ I xi'{x)dx +[{y/2) sup v{u)]'^ , 

ue(y/2,y) Jo «G(y/2,y) 

which is finite because |a;|i'(a;)c?x < oo (being X a bounded variation process). To obtain the bound on the first 

integral (jlip . we used the fact that the range for x is from to y/2 and, hence, the maximal range for u in the inner 
integral is u € [y/2,y]. Similarly, by Fubini's theorem, 

oo py^x pO ry—u 

v{u)v{x)dudx = / / v{u)v{x)dxdu 

) J — oo J y 

/O ry~u p-l ry-u 

I v{x)dxv{u)du + / v{x)dxv{u)du 

-1 J y J —oo J y 

/O p — 1 poo 

{—u)v{u)du+ / v{u)du I v{x)dx<oo. 
-1 J -oo Jy 

In the following proposition, we use Theorem 12.11 to establish a small-time estimate for the price of an out-of-the- 
money call option under the model in ([7|: 

Proposition 2.2 Assume that 



(i) / e^v{x)dx < oo and (ii) sup e^vlx) < oo, (12) 

J\x\>l \x\>e 

for any e > 0. Then we have the following small-time expansion for the price of a call option with strike K > Sq 
1 f°° 1 

-E{St-K)+^So {e'' -e'')+iy{x)dx + -So[d*2{k)-e''d2{k)]t + o{t) {t ^ 0) , (13) 



where k ~ log > is the log-moneyness and d2{k) = d2{k] b* , a, v*) with b* and v* given by 

v*[x) = e'^v{x) and b* = b + / x {e'' - 1) :y{x)dx + a'^ . (14) 

J\x\<l 

Remark 2.2 This result sharpens the asymptotic behavior ©, established by Levendorskii [LevOSj for a class of 
multi-factor Levy and Markov models under certain technical conditions. As explained in the introduction, for a Levy 
model, these conditions were relaxed by Roper [RoplO| and Tankov |TnkvlO| . Note also that, by imposing that i' has 
a positive Levy density, we are precluding the Black-Scholes case where there is a non-zero diffusion component with 
volatility cr and zero jump component, for which the implied volatility is just constant and equal to a. 

Proof. Without loss of generality, we assume that (Xt) is the canonical process Xt{uj) — uj{t) defined on = 
D([0, oo),]R) (the space of right-continuous functions with left limit lu : [0, oo) — >■ M) and equipped with the u— field 
J- = a{Xs '■ s > 0) and the right-continuous filtration J^t '■= ns>t(T(X„ : u < s). Following the density transformation 
construction of Sato jSat99] (see Definition 33.4 and Example 33.4 therein) and using the martingale condition we 
define P* on (51, J") such that 

P*(B) =E(e^nB), (15) 
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for any t > and _B G J-j. As explained in the introduction, we can interpret P* as the martingale measure associated 
with using the stock price as the numeraire. 

Let us first note that the price of a call option can be decomposed as follows: 

E{St~K)+ = E{Stls,>K)~KF{St>K) 

= 5'oE(e^' 1s,>k) - Soe^FiXt > k) 

= 5oP*(Xt > fc)-S'oe'=P(Xt > fc) (16) 

One can check that [Xt) is a Levy process under P* with characteristic triplet (6*, ct^, v*). For this result see the more 
general Theorem 33.1 in SatolSat99J. Finally, applying Theorem 12.11 to the probabilities under P and P* in (fT6|). we 
have 

- E(5t - K)+ = So f e''v{x)dx -K f v{x)dx + l-Sod*{k)t - ]-Kd2{k)t + o{t) ^ 0) , (17) 

* Jk Jk 2 2 

which simplifies to ■ 

Remark 2.3 Let us note that for a bounded variation process, the drift of X under the Share measure P* is the same 
as the drift under the measure P. Indeed, denoting 6q the drift under P* . we have that 



5o = &* — / xi'*{x)dx = b+ / X {e^ — 1) i'{x)dx — / xe^v{x)dx = h^. 

J{\x\<i} J\x\<i h\==\<^} 

Also, note that the call price approximation (|13p is independent of h. Indeed, let 

R{y- v) := d2(fc; 6, a, v) - {-a^iy'iy) + 2bv{y)) , 

which depends only on v as seen from the expression of ^2 in ([9]) . Then, using ([T4| , the second order term in ([T3| can 
be simplified as follows 

ai(fc) := ai{k]h,<7,v) := ^[d2{k;b* ,a,i^*) - e''d2{k;b,a,i^)] 

2 f -1 

= ^e''iy{k) + e'^vik) / a; (e^ - 1) v{x)dx + - [i?(fc; v*) - e''R{k; ly)] , 
which does not depend on b. The previous expression also shows that 

^2 

ai{k; b, a, v) — ai(fc; &, 0, v) = — e'"V(fc), 
and, hence, a non-zero volatility of a has the effect of increasing the call price approximation by ^-^e^v{k). 

2.1 Implied volatility 

Let (Jt{k) denote the Black-Scholes implied volatility at log-moneyness k and maturity t with zero interest rates, and 
let V{t,k) = at{k)'^t denote the dimensionless implied variance. Let 



poo 1 

(fc) := / (e^ - e'=)+i^(dx) and ai{k) := -[d;{k) - e''d2{k)] (18) 

J -OD 2 



denote the (normalized) leading order and correction terms in ()13p . By put-call parity, the dominated convergence 
theorem, and the stochastic continuity of the Levy process (Xt), we have 

limEiSt - K)+ = {So - K)+ , 

and from this we can show that V{t, fc) — > as t — t' 0. The following corollary shows more precisely how V{t, fc) — > 
as t and, hence, sharpening a result in Tankov [TnkvlOj (Proposition 4 therein): 
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Theorem 2.3 For the exponential Levy model in we have the following small-time behavior for the implied variance 
V{t, k) fork>0 

V{t,k)^Voit,k)[l + Viit,k) + oi-\)] (i^O), 



where 



Vo{t,k) = 



Vi{t,k) 



lu2 
2 

1^ 



log(i) 



■ los 



40rao(fc)' 



log ( - 



3/2' 



Proof. Sec Appendix \K\ ■ 

Remark 2.4 Multiplying (fT9|) by we have the following expansion for the implied volatility 



tlog(i) 



(19) 



(20) 



(21) 



and we see that (fc) — j- od as i — > 0+, as is well documented in e.g. CarrfcWu }CW03j (see also Roper [R"oplO| and 
|TnkvlO] ). The leading order term agrees with that obtained in Tankov [TnkvlO| and, moreover, we see that 

[t\og{^)]hat{k)^\k\/V2, (i->o), 

so the (re-scaled) leading order implied volatility smile is V-shapcd and independent of v, except that we require v to 
be non-zero. 



Remark 2.5 V{t,k) = 0{ . 1 ), so V{t,k) — > but slowly; in fact slower than fP for any p > 0. In particular, for a 

log I 

given desired "precision" bound e ^ 1, we will need t ~ 0{e~^/^) to ensure that V{t, k) — 0(e) and for the t-^-t- error 

term in ([T^ to be 0(e). For this reason, the call option estimate is more useful than the implied volatility estimate 
(|2ip in practice. We remark that in Corollary 8.3 of the very recent article by GaofcLee [GLllj . the authors give an 
expansion which sharpens (|19p . but proving their result is more involved and requires several preliminary lemmas 

Remark 2.6 Based on high-frequency statistical methods for Ito semimartingales, several empirical studies have 
statistically rejected the null hypothesis of either a purely-jump or a purely-continuous model (see, e.g., |AJ09bj . 
|AJ10j . |BNS06| ). If this really is the case, then our results show that theoretically, the small- maturity smile must 
tend to infinity, if put/call options are priced correctly. Nevertheless, this effect is often obscured in reality by market 
practicalities - high bid/offer spreads, daycount/settlement conventions, and times when the market is closed. However, 
even if we cannot trade an option with infinitesimally small maturity in practice, we can still look at rate at which 
the implied volatility smile steepens as the maturity goes small; typically it is difficult to fit the one of the fashionable 
class of purely continuous models (e.g. Heston, SABR, and other local-stochastic volatility hybrid models) to this kind 
of data, with realistic parameters. CarrfcWu }CW03j 's study of S&P 500 option price data (in contrast to the previous 
statistical approaches) also suggests that the sample path of the index contains both continuous and discontinuous 
martingale components (working under a risk neutral measure), and that, while the presence of the jump component 
varies strongly over time, the continuous component is omnipresent. 

In the same vein, Ait-SahaliafcJacod [AJ09a| define a jump activity index to test for the presence of jumps, which 
for a Levy process coincides with the Blumenthal-Getoor index of the process. |AJ09aj also proposes estimators of this 
index for a discretely sampled process and derive the estimators' properties. These estimators are applicable despite 
the presence of a Brownian component in the process, which makes it more challenging to infer the characteristics of 
the small, infinite activity jumps. When the method is applied to high frequency stock returns, jAJ09aj found evidence 
of infinitely active jumps in the data and they were able to estimate the index of activity. 



7 



3 Time-changed Levy processes 



3.1 A formula for out-of-the- money call option prices 

In addition to the Levy process (Xt) of Section [2l we now consider a random clock (Tj) defined on {^,¥,T) and 
independent of X. A random clock is a right-continuous non-decreasing process such that Tq = 0. We consider a 
time-changed Levy model of the form 

5t:=5oe^S with Zt := Xt^. (22) 
As explained in the introduction, this type of model is important because it can incorporate volatility clustering effects. 

Given that e-^* is a martingale under P (relative to the natural filtration generated by X), it is known that (^f) 
above is a martingale under P relative to the natural filtration generated by the random clock Tt and the time-changed 
process Zt (sec Lemma 15.2 in |CT04| ). Note also that our simplifying assumption ([8|) implies that 

^' = ^"E(^ ^^^^ 

because E(e^*) — 1. jCGMY0"3] (Section 4.2) shows that the price process ((23)) is free of static arbitrage opportunities. 
Furthermore, under certain conditions (e.g. if X has infinite jump activity and (Tt) is continuous), (t(T„ : u < t) G 
<j{Xt^ : u <t) (see, e.g.. Theorem 1 in |Win01j ). and hence ([22]) will be a martingale relative to the filtration generated 
by only the time-changed process {Zt) or, equivalently, the filtration generated by the stock-price process {St). In that 
case, the model will be free of dynamic arbitrage opportunities by the sufficiency part of the First Fundamental 
Theorem of Asset Pricing. 

Let N be the set of P-nuU sets of F and define a probability measure P on := a{Zt,Tt : t > 0)\/ M such that, for 
any t > 0, _ 

P(B) =E(e^'lB) , (24) 

whenever B ^ Ft a{Zu,Tu ■ u < t) V Af. We note that P is well defined since {e^^}t>o is a P-martingalc relative to 
{Ft}t>o- The following proposition will play a key role in the sequel: 

Proposition 3.1 Suppose that the assumptions of Proposition are satisfied and let {b* ,a'^,i'*) be defined as in 
Then, under P, the process {Zt) in h22\l has the same distribution as a Levy process with the characteristic triplet 
{b*,a^,i^*) evaluated at the independent random clock Tt. 

Proof. Fix = tQ<---<tn— t<oo and Mi, . . . , u„ G M. Then, using the independence between T and X, 

n n n 

E(exp{*^u,(Zt^ -Zt^_J}) =E(exp{Zt + i^u,(Zt^. -Zt^_J}) =E(exp{^i(7., -z)(Xt,^. -Xt,^_J}) 

n n 

= E(exp{^(Tt^. - Tt^_,)^P{u, - z)}) = E(exp{^(Tt^. - Tt^_jr{u,)}) ■ 
i=i j=i 

The last expression corresponds to the characteristic function of a process of the form X^^, where {Xt) is a Levy 
process with triplet {b* ,a'^, v*) defined on (ri,P, J") and independent of the random clock {Tt). ■ 
In light of the previous result, we have the following representation for call option prices: 

¥.{St-K)+ = E{Stlst>K)-KnSt>K) 

= SoE{exp{Zt)lzt>k) - Soe''¥{St > K). 

= S^¥{Zt>k)-Soe''nZt>k). (25) 

We emphasize again that, under P, Zt has the same distribution as a Levy process with characteristic triplet {b* ,a'^,i/*) 
evaluated at an independent random clock {Tt). Hence, as for the pure-Levy model case, the problem of finding 
small-time expansions for out-the-money option prices reduces to finding small-time asymptotics of the corresponding 
distribution functions. 
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3.2 Small-time asymptotics for the time-changed Levy model 

In this section, we determine the asymptotic behavior of out-the-money call option prices. We consider random clocks 
{Tf) that are absolutely continuous with non-negative rate process (Yt) (i.e. Tt = /p Y^ds) such that Yq > 0. We will 
also refer to the following conditions in the sequel: 

(i) EYt-EYo = 0{t), (ii) limsup EF^^ < oo, (iii) \im - [EYt - EYq] ^ e [0 , oo) , (26) 

(iv) limsup EY,^ < oo, (v) lim \eT^ ^ pe(0,oo). (27) 
t\,o tN.0 

In the case that (Yt) is a stationary process with finite moment of third order, KY^^ is constant for k = 1,...,3 
and (i)-(iv) are automatically satisfied. Also, if Yt Yq and (iv) are satisfied, then (v) holds true with p = EYq. 
Indeed, note first that liuis^o^^s = since (F(^)t<to are uniformly integrable for small enough to by (iv) above. 
Also, since /t^ < Y^ds/t (by Jensen's inequality) and limf_^o E /q Y^ds/t = Elimt^o /q Y^ds/t, so the dominated 
convergence theorem implies that 

lim ^ET? = E lim (- [ Y^ds] = EY^ . 

The following result gives the small-time asymptotic behavior of the tail distributions of time-changed Levy models: 
Theorem 3.2 Suppose that the conditions of Theorem \2.1\ are satisfied as well as conditions (i)-(ii) of i86\) . Then, 

ViZt > x) = tEYoiy[x, oo) [1 + 0{t)] , {t 0). (28) 

If additionally, conditions (iii)-(v) of |g7| j are satisfied, then 

V{Zt > x) = tEYQi/lx, oo ) + 7;{pd2{x)+ji^[x,oo))t^ + o{t^), (i->0), (29) 

where d2 is the same as in Theorem \2.1[ 
Proof. See Appendix [X] ■ 

Remark 3.1 A very popular rate process in applications is the Cox- IngersoU- Ross (CIR) diffusion process, defined by 

dYt^ K{e ~Yt)dt + ay/YtdWt, (30) 

where (Wt) is a standard Brownian motion, Yq is an integrable positive random variable independent of W, and 
K,9,a > are such that K9/a^ > 1/2 (which ensures that y = is an inaccessible boundary). If Yq ~ r{^^, ^), 
the proces (Yt) is stationary and EY^'' is finite and constant in t for any fc > 1. In particular, (i)-(v) arc satisfied with 
p = EYq. In the non-stationary case, it is known that EYt — EYq = {6 — EYq) (1 — e^"*) and (i) & (iii) are satisfied 
with 7 = k{9 — EYq). The other conditions in (|26ll27p will also hold true. Thus we conclude that the time-changed 
Levy model with CIR speed process satisfies: 

P{Zt > x) ^ tEYQiylx, oo) + (EYq d2{x) + n {e~EYQ)v[x,^))-t' + o{i'), (t^O). 

We are now ready to give the small-time asymptotic behavior of out-the-money call option prices and the corresponding 
implied volatility: 

Corollary 3.3 Under the conditions of Proposition [37^ we have the following small-time expansions 

lE{St-K)+ = SQEYoaQik) + 5o [pai(fc) + 7ao(fc)] ^ + o(i) (< ^ 0) , (31) 
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where k ~ logK/So > and qq, ai are the first and second order terms of the call price approximation il3\} as defined 
in 118\) . Furthermore, we have the following small-time behaviour for the implied variance V{t, k) for k > 



Vit,k) = Voit,k)[l + Viit,k) + oi-\)] 



where 



Voit,k) 
Viit.k) 



"2 ^ 

1^ 



log(i) 



lo{ 



/ 4v^E(yo)ao(fc)e-''/^ 



(t^O), 



(32) 



log - 



1 3/2^ 



(33) 



Proof. The expansion ((3T|) follows from the representation (p5)) and ((29| . The asymptotics (|32p follows from the proof 
of Thorcm [ 



Remark 3.2 As it was indicated before, the time-changed Levy model ([22l) was introduced to account for the volatility 
clustering exhibited by financial time series. Indeed, the process {Yt)t controls the speed of the random clock so that 
when Yt is high, the random clock runs faster and, hence, the price process exhibits more variability. Another approach 
to incorporate stochastic volatility is via stochastic integration along the lines of the following jump-diffusion model 

d\n{St/SQ)=fi{Yt)dt + a{Yt)dWl+dZt, dYt = a{Yt)dt + -f{Yt)dWt^^\ (34) 

where M^*-^-* and M^*-^^ are two (possibly correlated) Brownian motions and Z is a pure-jump process. For a comparison 
of these two methods, we refer the reader to Chapter 15 of |CT04| . Recently, |FLGHll| have provided small-time 
expansions for vanilla option prices under the stochastic model p4p when Z is a pure-jump Levy process independent 

of y. 



4 Small-time, small log-moneyness asymptotics 

In this section, we survey the behavior oiF{Xt > k) for a Levy process AT, when t — and k = kt also converges to zero 
at an appropriate rate. We can think of this scaling as a small-time, small log-moneyness regime. As an application, 
we deduce the asymptotic behavior of at-the-money call option prices for a CGMY model. 

4.1 Levy models with non-zero Brownian component 

Several financial models in the literature consist of a Levy model with non-zero Brownian component. The most 
popular models of this kind are the Merton model and Kou model determined by the characteristic functions 

E(cxp(mAi)) cxp[t{ibu - -a^u'^ + iuX {- — — —))], 

2 A-i- — TO A_ + TO 

E(cxp(TOAt)) = cxp[t{ibu - icr^u^ + A (g-'^'^Va+v" _ i))]. 

It turns out that, for a general Levy process (Xt) with cr 7^ 0, 

liniE(exp(TO,At/%/t)) = cxp(-i(T^u^) , 

(see e.g. pp. 40 in |Sat99| for a formal proof). The right-hand side is the characteristic function of a Normal A^(0, cr^) 
random variable Z, thus [Xt/y/t) converges weakly to a Normal distribution with variance cr^ and 

limP(At/Vi > x)= ¥{Z > x) . 
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4.2 The CGMY model and other tempered stable models 



The so-called CGMY model is a pure-jump Levy process determined by a Levy density of the form 

= ^i+Y ^{^^>o} + ^{^<o}- (35) 

for C, G,M > and Y € (0,2). As explained in the introduction, the CGMY model is a particular case of the more 
general KoBoL class of models, named after the authors |Kop95| (who first introduced the symmetric version of the 
model under the name of "truncated Levy flights") and |BL02| . The term CGMY was introduced later on by Carr et al. 
|CGMY0"2| . This process is a tempered stable process (see Section 4.5 in ContfcTankov [CT04| ) . and its characteristic 
function is given as 

0t(u) = E(e™^*) = cxp \tCT{-Y) {{M - iuf + (G + iuf - M^' ~G^'}+ ibutl , (36) 

for y ^ 1 and some constant & G M (see [CT04| for the formula when y = 1). We note that we must have M > 1 for 
(|12p to be satisfied, and under this condition, X is again a CGMY process under P* with parameters C* — C, Y* ^ Y, 
M* = M — 1, and G* = G -I- 1. In the bounded variation case {Y < 1), b coincides with the drift bo. 

The following result characterizes the small-time behavior of V{Xt > kt) with small log-moneyness kt ^ xt^^^ . 

Proposition 4.1 For the CGMY model with Y £ (1,2), (Xt/t^/^) converges weakly to a symmetric Y-stable distri- 
bution as i — > 0. Concretely, 

liniP(Xt/t^/^ > a-) = ¥{Z > x) , 

where Z is a symmetric Y-stable random variable with scale parameter c ~ (2Gr(— y)| cos(^y7r)|)"'^/^; i.e. Z has 
characteristic function 

Ciu) = exp(-2Gr(-r)| cos(ir^)| \u\^) . 

Remark 4.1 Note that Z has infinite variance because Y < 2. The stable distribution was famously used by 
Mandelbrot jMan63| to model power-like tails and self-similar behaviour in cotton price returns. 

Proof. Let 

ip{u) = Gr(-y)((M - iu)^ + (G + iuf - - G^) + iub (37) 
denote the characteristic exponent for the CGMY process. Then we have 

C(u) = limexp(tV(^)) = exp(-Gr(-y)|(-z)^ + I^D , 

where we used that Y S (1,2). Q(u) is continuous at zero and we recognize C,(u) as the characteristic function of a 
symmetric alpha-stable distribution. Thus, by Levy's convergence theorem (see Theorem 18.1 in Williams jWill9 1 j ) . 
the sequence of random variables [Xt/t^/^) converges weakly to Z. The second result follows from the Lemma on 
page 181, chapter 17 in jWill91j . ■ 

Remark 4.2 Proposition 14. II is a particular case of a result shown in Rosihski jRos07| where a more general class of 
tempered Levy measures is considered. Concretely, |Ros07j considers Levy measures of the form 

y{A) ^ f f lA{uw)u-^-^e-''duR(dw), (38) 
Jr Jo 

for a measure R such that i?({0}) = and J^i\w\^ A \w\^)R{dw) < oo. The CGMY model is recovered by taking 
R{dw) = CM^ 5{M-^}{dw) -i- CG^ 6{_Q-i^{dw). In light of Rosinski's Theorem 3.1, it follows that Proposition |4TT] also 
holds true for Y £ (0, 1) (finite- variation case) provided that {Xt) is driftless, i.e. b in must be (otherwise, we 
have to replace Xt by Xt — bt) . Note that under P* , X is also driftless (see Remark 12. 3p . 
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Another well-know class of Levy processes is the Normal Inverse Gaussian (NIG) model, introduced in Barndorff- 
Nielsen }Bar97| . for which the characteristic function is given by 

E(exp(mXt)) = exp[^t6{^/^!^~~(f+iu)^ y/a^ - 13^)] . 

The Levy density of the NIG model takes the form j/(.t) = Ce^'^ Ki{B\x\) /\x\ where Ki is the modified Bessel function 
of second kind and A,B, and C are certain positive constants (see jCT04j for their expressions). Hence, one can view 
the NIG process as an improper tempered stable process in the sense of Rosiriski jRosOTj . It is also easy to see that 

lim E(exp(mXt/t)) ~ cxp[—tS\u\] . 

The right-hand side is the characteristic function of a symmetric alpha-stable random variable Z with a = 1 and 
scale parameter S i.e. a Cauchy distribution; thus by the same argument we see that (Xt/t^^^) converges weakly to a 
symmetric Cauchy distribution: 

\iTjiF( Xt/t^/^ >x)= F(Z > x) . 

t-^0 

4.3 At-the-money call option prices for the CGMY model 

Our approach to deal with at-the-money call option prices is based on the following result from CarrfcMadan [CM09j : 

^EiSt-K)+=¥*iXt-E>\og^), (39) 

where E is an independent exponential random variable under P* with parameter 1. Now set K = Sq- Consider the 
CGMY model with Y € (1, 2). The idea is to use the small-time, small log-moneyness result in the previous section. 
Indeed, note that 

/>oo /-oo 

t-yY^*{Xt>E)=t-^/^ e-^F*{Xt>x)dx^ e-''^''^F*{Xt>t^/^u)du. (40) 

Jo Jo 

From our Proposition 14. 1[ 

P*{Xt > t^l^u) ¥*{Z > u), 
for any u > 0, where Z is a symmetric a-stable r.v. under P*. The previous fact suggests the following result: 

Proposition 4.2 Suppose that X is a CGMY process under ¥ with Y G (1,2). Then, the at-the-money call option 
price has the following asymptotic behavior: 

limt-^^^E{St - So)+ = SoE*{Z+), (41) 
where Z is a symmetric Y -stable r.v. as in Proposition 4^.1. 

Proof. See Appendix 1X1 ■ 

In order to justify the previous argument, we will need the following estimate: 

Lemma 4.3 Let X denote a symmetric CGMY process under V (hence G ~ M) with Y <E (1,2), M > 1, and C > 0. 
Then, there exists a universal constant K > such that 

P* {Xt >x)< Kx-^t. (42) 

for any t > and a; > satisfying t{b + /|2|<2./4 ^{^^ ~ < x/4. 

Proof. See Appendix El ■ 
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Remark 4.3 As seen in the proof of Lemma l473l the estimate (|42p is vahd for any pure-jump Levy process admitting 
a symmetric Levy density i'{x) such that 

^-M\x\ 

for some Y e (1, 2), C > 0, and A/ > 1 . Moreover, as seen in the proofs of Proposition 14.21 if we further assume that 

(t-^/^^t)^ A iZt\, (43) 

as i ^ under P* (for a symmetric F-stable process {Zt)t), then the asymptotic behavior (|4ip wih also hold. Condition 
(|43p holds for a wide range of processes (see, for instance, Proposition 1 in jRTllj for relatively mild conditions). 



4.4 At-the-money implied volatility 

Proposition 4.4 For the CGMY model with y G (1, 2) m Proposition \4-. S\ we have the following small-time behaviour 
for the at-the-money implied volatility at(0) 

lim<Tt(0)/t^/'^"^ = V2^E*{Z+) . 

Proof. We first recall that the dimensionless implied variance V{t, 0) — d-t{0)^t — > as t 0. Equating prices under 
the the Levy model and the Black-Scholes model, we know that for any S > 0, there exists a. t* = t*{S) such that for 
all t <t* we have 



E*(Z+)ti/^(l-^) < ^E{St-So)+ < ^^^^\ l + S). 

Do v Stt 



Re-arranging, we see that 



1-s ^ VW^ 



l + S - V2^E* {Z+)t^/^ ■ 
We proceed similarly for the upper bound. ■ 

5 Robust pricing of variance call options at small maturities 

Let {Xt) denote the general Levy process defined in section[21 The quadratic variation process [X]t = f ^i-|-X].s<t(^^s)^ 
is a subordinator and has Levy density given by 

^^^^ 2^ 2^ ' 

(see e.g. |CGMY0"5] ) . The function f{y) = [y — K)^ for K > satisfies the conditions of Theorem 1.1 in Figueroa- 
L6pez [FL08j . so we have 



oo 



\ E{[X]t - K)+ = I {y- K)+q{y)dy + 0{t) (t ^ 0) (44) 







oo 







iy-Kh[^ + ^^]dy + Oit) 



oo 

2 



{x^ - K)+iy{x)dx + 0{t) (45) 



— OO 



jE{Xf~K)+ + 0{t) 
iE[(ln£)2-X]+ + 0(i) {t^O). 
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From this we sec that an out-of-thc-moncy variance caU option of strike K which pays {[X]t — at time t is worth 
the same as a European-style contract paying ((In -1^)^ — K)j^ at time t as t — >■ 0, irrespective of y{-)- Note that the 
diffusion component of Xt does not show up at leading order for small t. We also remark that the higher order terms 
in and (P5)) can be obtained by using the expansions in Theorem 12.11 and the following identities: 



6 Numerical examples 

In their seminal work, Carr et al. }CGMY02| calibrated the CGMY model and the Variance Gamma (VG) model to 
option closing prices of several stocks and indices. In this section, we shall use some of their calibrated parameters to 
illustrate the approximation proposed in this paper. As in Section [2l we are assuming below that the risk-free rate r 
and the dividend rate q are both set to be zero. 

Using IBM closing option prices on February 10th, 1999 and maturities of 1 and 2 months, |CGMY0"2] report the 
following calibrated parameters for the VG model: 



where ct, v, and are the three parameters characterizing the VG process (see e.g. |CT04| ). and 77 is the volatility of 
an additional independent Wiener component. In order to assess the accuracy of the call price approximation ()13p . 
we have plotted (in Figure [ij the first and second order approximations of E{St — K)+/t as a function of the log 
moneyness k = \ogK/S[) for S'o = 1 and time-to- maturities t = 5/252 and t = 10/252 (in years). We have also plotted 
the "true" option prices obtained via an inverse Fourier Transform (IFT) method (see Theorem 5.1 in |Lee04| for the 
case G = Gi corresponding to the call option payoff with a > 0). Table [T] also shows the numerical approximations for 
1000 X E(S't — K)+/t corresponding to four maturities, together with the numerical values obtained via the IFT. Note 
that the first order approximation (i.e. 1000 x /_^(e^ — e^)+v{x)dx) is independent of time-to-maturity t. The graphs 
show that the second order approximation significantly outperforms the first order approximation. The corresponding 
table shows that the second order approximation is quite good for maturities of 5 to 10 days and logmoneyness values 
larger than 0.1. 

The numerical values via the IFT method were implemented in Mathematica. while the coefficient Q was computed 
using numerical integration routines of Mathematica. This computation is typically slow due to the singularity of the 
Levy density v and the cumbersome double integrals. A much faster numerical method, valid for bounded variation 
Levy processes, is described in [FLlOj (see below for an illustration of this method). 

In order to illustrate the performance of the approximations for larger volatility values, we now consider the 



which were calibrated to fit INTEL option data as reported in |CGMY0"2] . The results are shown in Figure [T] for 
S'o = 1 and time-to-maturities t ~ 5/252 and t = 10/252 (in years). Table [2] shows the numerical approximations for 
1000 X E(S't — K)j^/t corresponding to four maturities. We also show the numerical values obtained via the IFT. The 
second order approximation is again quite good for mid-range log-moncyness values and no noticeable difference is 
observed even though 7/ is significantly larger. 

For the case of Microsoft option prices on December 9th, 1999 and maturities of 1 and 2 months, |CGMY0"2| report 
the following parameters for a CGMY model: 



Table [3] shows the numerical approximations for 1000 x E(S't — K)+/t corresponding to four maturities, together with 
the numerical values obtained via the IFT (computed using Mathematica). As before, the approximations perform 
quite well and we are able to attain a decent approximation even for a maturity of 20 days. To compute the second 
order approximations (or more specifically, to compute the coefficient (O), we have employed the method in |FL10| . 




cr = 0.4344, V = 0.1083, 6 



.3726, 77 = 0.0051, 



parameters: 



(7 = 0.1452, 61 = -0.1497, 7^ = 0.1536, 77 = 0.0869, 



C=l.l, G = 5.09, A/ = 8.6, F = 0.4456. 
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Figure 1: Here we have plotted the leading order term (grey line) and the correction term (solid blue line) of the 
approximation (fT3|) for jK{St — K)+ as a function of the log-moneyness x = k = log K/Sq for a Variance Gamma 
model with an independent Brownian component. The parameters of the VG model arc a = 0.4344, i/ = 0.1083, 
and = —.3726, while the volatility of the independent continuous component is 77 = 0.0051. Left and right panels 
corresponds to the expiration times t = 5/252 and t = 10/252, respectively. The numerical "true" option prices 
obtained via the IFT are also shown (dashed grey line). 
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X 


1st 
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0.05 


234.6977 
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258.4404 


254.5295 


282.1831 
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329.6684 


277.3445 


0.06 
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244.4061 


0.07 
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148.9089 
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116.7270 
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62.1087 
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74.3137 


73.9493 


86.5186 


84.8855 


110.9285 


103.4274 


0.14 


53.3465 


55.5346 


55.5269 


64.2872 


64.0541 


75.2279 


74.1246 


97.1093 


91.6844 


0.15 


45.9096 


47.8674 


47.8636 


55.6984 


55.5669 


65.4873 


64.7996 


85.0649 


81.3079 


0.16 


39.5787 


41.3278 


41.3269 


48.3238 


48.2701 


57.0689 


56.7045 


74.5590 


72.1326 


0.17 


34.1752 


35.7358 


35.7372 


41.9783 


41.9835 


49.7815 


49.6660 


65.3878 


64.0145 
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29.5521 


30.9433 


30.9463 


36.5080 


36.5571 


43.4639 


43.5376 


57.3758 
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25.5884 


26.8275 


26.8317 


31.7841 
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38.1947 
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50.4628 


0.2 


22.1834 


23.2864 


23.2913 


27.6985 


27.8019 


33.2136 


33.5313 


44.2438 


44.8227 



Tabic 1: Approximations (US]) for 1000 x iE(S't - K) + as a function of the log-moneyness x = k = log K/ Sq for a 
Variance Gamma model with an independent Brownian component. The parameters of the VG model are a = 0.4344, 
1^ = 0.1083, and 9 = —.3726, while the volatility of the continuous component is 77 = 0.0051. The column "1st" indicates 
the first order approximation (which is independent of t) . The column "2nd" refers to the second order approximation 
term. 
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Figure 2: Here we have plotted the leading order term (grey line) and the correction term (solid blue line) of the 
approximation (fT3|) for ■|-E(5't — if)+ as a function of the log-moneyness x = k = log K/ So for a Variance Gamma 
model with an independent Brownian component. The parameters of the VG model are a = 0.1452, 9 ~ —0.1497, 
v ~ 0.1536, while the volatility of the independent continuous component is 77 = 0.0869. Left and right panels 
corresponds to the expiration times t = 5/252 and t = 10/252, respectively. The numerical "true" option prices 
obtained via the IFT are also shown (dashed grey line). 
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1.0608 


1.1027 


1.3485 


1.5239 


1.9241 


2.6763 


0.13 


0.5298 


0.5698 


0.5709 


0.7297 


0.7598 


0.9297 


1.0562 


1.3295 


1.8801 


0.14 


0.3643 


0.3922 


0.3930 


0.5037 


0.5252 


0.6430 


0.7343 


0.9216 


1.3240 


0.15 


0.2513 


0.2708 


0.2714 


0.3486 


0.3641 


0.4460 


0.5119 


0.6406 


0.9344 


0.16 


0.1738 


0.1875 


0.1879 


0.2420 


0.2531 


0.3101 


0.3577 


0.4464 


0.6607 


0.17 


0.1205 


0.1301 


0.1304 


0.1683 


0.1763 


0.2161 


0.2504 


0.3117 


0.4679 


0.18 


0.0837 


0.0905 


0.0907 


0.1173 


0.1231 


0.1509 


0.1757 


0.2181 


0.3318 


0.19 


0.0583 


0.0630 


0.0632 


0.0819 


0.0861 


0.1056 


0.1234 


0.1528 


0.2356 


0.2 


0.0407 


0.0440 


0.0441 


0.0573 


0.0603 


0.0740 


0.0869 


0.1073 


0.1675 



Tabic 2: Approximations (US]) for 1000 x iE(S't - K) + as a function of the log-moneyness x = k = log K/ So for a 
Variance Gamma model with an independent Brownian component. The parameters of the VG model are <j = 0.1452, 
= —0.1497, V — 0.1536, while the volatility of the continuous component is 7/ = 0.0869. The column "1st" indicates 
the first order approximation (which is independent of t) . The column "2nd" refers to the second order approximation 
term. 
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Timc-to-mat. t 


1/252 


5/252 


10/252 


20/252 


X 


1st 


2nd 


IFT 


2nd 


IFT 


2nd 


IFT 


2nd 


IFT 


0.05 


118.8662 


120.2883 


120.5386 


125.9768 


125.9179 


133.0875 


131.5844 


147.3088 


139.5891 


0.06 


99.6004 


100.8808 


101.1351 


106.0023 


106.0868 


112.4042 


111.5177 


125.2081 


119.9024 


0.07 


84.3149 


85.4610 


85.7023 


90.0455 


90.1924 


95.7760 


95.2726 


107.2372 


103.5827 


0.08 


71.9095 


72.9321 


73.1727 


77.0226 


77.2339 


82.1358 


81.9201 


92.3620 


89.9114 


0.09 


61.7191 


62.6303 


62.8747 


66.2750 


66.5275 


70.8309 


70.8150 


79.9426 


78.3608 


0.1 


53.2682 


54.0799 


54.3141 


57.3264 


57.5892 


61.3846 


61.4910 


69.5011 


68.5328 


0.11 


46.1664 


46.8892 


47.1192 


49.7805 


50.0626 


53.3947 


53.6011 


60.6229 


60.1205 


0.12 


40.1763 


40.8204 


41.0433 


43.3967 


43.6782 


46.6171 


46.8806 


53.0579 


52.8833 


0.13 


35.0705 


35.6445 


35.8690 


37.9408 


38.2302 


40.8111 


41.1241 


46.5517 


46.6292 


0.14 


30.7034 


31.2154 


31.4361 


33.2632 


33.5566 


35.8230 


36.1693 


40.9425 


41.2037 


0.15 


26.9570 


27.4140 


27.6311 


29.2418 


29.5285 


31.5266 


31.8864 


36.0962 


36.4806 


0.16 


23.7163 


24.1244 


24.3391 


25.7565 


26.0433 


27.7968 


28.1703 


31.8772 


32.3565 


0.17 


20.9085 


21.2731 


21.4858 


22.7315 


23.0167 


24.5545 


24.9355 


28.2005 


28.7454 


0.18 


18.4722 


18.7982 


19.0082 


20.1025 


20.3798 


21.7327 


22.1107 


24.9933 


25.5756 


0.19 


16.3432 


16.6349 


16.8407 


17.8017 


18.0761 


19.2602 


19.6377 


22.1771 


22.7868 


0.2 


14.4852 


14.7463 


14.9482 


15.7910 


16.0580 


17.0968 


17.4672 


19.7084 


20.3280 


0.21 


12.8531 


13.0870 


13.2891 


14.0226 


14.2859 


15.1920 


15.5580 


17.5310 


18.1563 


0.22 


11.4193 


11.6289 


11.8268 


12.4672 


12.7267 


13.5150 


13.8752 


15.6108 


16.2344 


0.23 


10.1595 


10.3474 


10.5434 


11.0990 


11.3517 


12.0385 


12.3891 


13.9176 


14.5312 


0.24 


9.0459 


9.2145 


9.4085 


9.8885 


10.1371 


10.7310 


11.0744 


12.4161 


13.0193 


0.25 


8.0621 


8.2133 


8.4040 


8.8179 


9.0625 


9.5737 


9.9096 


11.0853 


11.6753 


0.26 


7.1931 


7.3287 


7.4365 


7.8714 


8.1099 


8.5498 


8.8759 


9.9065 


10.4792 


0.27 


6.4212 


6.5430 


6.7291 


7.0301 


7.2645 


7.6389 


7.9573 


8.8567 


9.4132 


0.28 


5.7374 


5.8468 


5.8054 


6.2842 


6.5132 


6.8309 


7.1400 


7.9243 


8.4622 


0.29 


5.1285 


5.2267 


5.4878 


5.6194 


5.8445 


6.1103 


6.4118 


7.0920 


7.6128 


0.3 


4.5867 


4.6749 


4.8038 


5.0275 


5.2487 


5.4683 


5.7624 


6.3499 


6.8534 


0.31 


4.1050 


4.1842 


3.4559 


4.5009 


4.7173 


4.8968 


5.1826 


5.6886 


6.1739 


0.32 


3.6746 


3.7457 


3.7292 


4.0301 


4.2427 


4.3856 


4.6643 


5.0966 


5.5652 


0.33 


3.2905 


3.3543 


3.6098 


3.6097 


3.8185 


3.9289 


4.2006 


4.5673 


5.0195 


0.34 


2.9479 


3.0053 


3.2470 


3.2346 


3.4391 


3.5212 


3.7855 


4.0944 


4.5299 


0.35 


2.6410 


2.6925 


2.8716 


2.8983 


3.0991 


3.1555 


3.4134 


3.6701 


4.0903 



Table 3: Approximations (|13p for 1000 x ^£(5'^ — K)^ as a function of the log-moneyness x ^ k ~ logK/So for the 
CGMY model with parameter values C = 1.1, G = 5.09, M = 8.6, and Y = 0.4456. The column "1st" indicates the 
first order approximation (which is independent of t). The column "2nd" refers to the second order approximation 
term. 
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We now proceed to illustrate the performanee of the implied volatility approximations described in Section 12.11 
Concretely, we analyze the relative error of the approximations 

y^p), yMMi+MM)^ ,46) 

Let us first analyze the Variance Gamma model with parameter values as above. The left panel of Figure [3] shows 
the relative errors {fft.i ~ ^t)/o't and ((Tt,2 ~ ^t)/^t as a function of time-to-maturity t for values of k ranging from 
0.1 to 0.3. Note that both at^i and at,2 consistently underestimate the true implied volatility. For k = 0.3, the first 
order approximation is actually quite good with a relative error of about —5% uniformly in t and it is only for very 
small values t (less than 3 days) when at^2 is better than at^i- However, for the other values of fc, at.2 significantly 
outperforms (Jt,i- For instance, for k = 0.2, the relative error of at^i ranges from —19% to —34% with a mean absolute 
error of 27.0%, while the relative error of at,2 rages from —4.4% to —23% with a mean absolute error of 14.2%. The left 
panel of Figure |4] compares the term structure of the approximated implied volatilities to the "true" implied volatilitjU. 
The right panel of Figure [3] shows the analog results for the CGMY with parameter values as above. The results are 
qualitatively similar to those of the Variance Gamma model. However, all the approximations seem to perform better 
in terms of error stability in time and accuracy. For k — 0.2, the relative error of at,i ranges from —12% to —20% 
with a mean absolute error of 18.6%, while the relative error of CTt_2 ranges from 0.83% to —13% with a mean absolute 
error of 9.25%. The right panel of Figure |4] compares the term structure of the approximated implied volatilities to 
the "true" implied volatilit}^. 




Figure 3: Relative errors of the implied volatility approximations for the VG and CGMY models as function of time 
to maturity using the two estimators at^i and at^2 in (|46|). 
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A Proofs 

Proof of Theorem 

We know that V{t, k) — >■ 0. Equating call prices in the small-time limit under the exponential Levy model (using 
Proposition 12. 2|) . and the Black-Scholes model with zero interest rates and implied variance V = V{t,k) (using e.g. 
Proposition 3.4 in [FJLlOj or Lemma 2.5 in |GHLOW09] ) we know that for any (5 > 0, there exists a. t* = t*{5) such 
that for all t <t* 

tao{k){l-5) < ^¥.{St-K)+ < e-i'^'(i~'5)/^(*'fc) . (A-1) 
So 



Re-arranging, we see that 



-V{t,k)\og[taoik)il-d)] > h\l-d)^ 



or 

F(t,fc).log(i) > h^l - S) + V{t,k)log{l - 5) + V{t,k)\ogao{k) . 

V{t, k) — 0, so this yields a lower bound for V{t, k). Using a similar argument for the corresponding upper bound, 
we establish the leading order asymptotic behaviour for the implied variance as 

1^,2 

V{t, k) - Vo{t, k) := {t^O). (A-2) 

Now let V{t,k) = Vo{t,k)[l + Vi{t,k)] and note that Vi{t,k) = o(l) as i — > 0. Then for any (5 > 0, there exists a 
t** = t**{S) such that for t < t** we have 

III ^E{St - K)+ - ao{k)] - ai (fc) > -S . (A-3) 

t t Do 

Re-arranging, we have 

too(fc) + (ai(fc) - S)t^ < ^E{St - K)+ . (A-4) 

So 

Using this bound and again equating small-time call prices under the Levy model and the Black-Scholes model, we 
have that there exists a positive constant c such that for t small enough 

tao{k) + {a,ik)-S)t' < ^E{S, -K)+ < £l!plMe-|'=Vv(*.fc) (1 + e^t, fc)) 

oo yzn k'' 

^ e^^yo(t,fc)i(l-fT/i(f,fc))i 

V2^fc2 ^ ^' " 



^ £i!Ml^e-5'='/{^«(*'^-)(i+^^(*''=»>(l+f(<,fc)), {t-,0). 
\/27r A:^ 
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where := (1 + Vl(i, A:))^/^(l+cy(i, fc)) — l, which converges to as t ^ 0. Dividing both sides by i = e 2'=VVo(t,fc) 

we have 

ao(fc) + (ai(fc) - S)t < g^fc^v,(t,fc)/{y„(^fc)(i+v-,(t,fc))}(^ ^ ^^^^ ^ 

and re-arranging we obtain 



Viit,k) 
1 + Vi{t,k) 



2 



> ^Vo{t, k) log [(ao(fc) + (ai(fc) - 5)t)V2^k^e---''Vo{t, fc)" V(l + ^(i, k))] 



St 



-Vo{t,k)\og [(ao(fc)+tal(fc))^/2^fcVsfcyo(t,fc)-i] +-yo(t,fc)log[ ^ 

ai(fc) 



log(i) 



■ los 



l^^^^^M^[iog(i)]i 



log(i 



los 



ao(fc) 



1 - 



ao(fc)+ai(fc)t 



l + f(i,fc) 



Vi(t,fc) 



Note that = fc) = O (^ '°|;'°|' ^ and £' = f '(t, fc) = o(j^) since fc) -> as t ^ 0. Solving the inequality 



Vi 
l + Vi 



> Vi we find that 



Vi+£' 



, Vi + 2^1^ + e 



12 



l-{Vi+£') l-{Vi+£') 

Since £'{t, k) > Vl{t, k) for t sufficiently small, we conclude that Vi > Vi+o{-r-^)- Proceeding similarly for the upper 
bound, we conclude that 



Vi{t,k) = Vi{t,k)+o{—^). 

log* 



as t — > 0. 



Proof of Theorem [mi Let F{t) := F{Xt > x) and B := iy[x, oo) + sup^^o F{Xt > x)/t. In the light of Theorem 
2.11 there exist constants to > Q and A' < oo such that 



"{Xt > x) - v[x,oo) 



< Kt, 



for any < t < to. Next, conditioning on Tt, 



jPiXT,>x)^jEFiTt)='^^^^ETt + jE 



Tt 



F{Tt) - v{x, c5o) \ Tt 



Let Rait) denote the second term on the right-hand side, which we can bound as follows: 



1 



m < - E( i{T,<to} 



F{Tt)-iy[x,^) 



1„ , „1 



^ ^2 I THi frfi2\ 



F{Tt)-iy[x,^) 



Tt 



< K-ETt + B-E{ l{T,>t,}Tt) < K-ETl + --E(r/) 

I I T to t 

using a Chebyshev upper bound. Combining the previous bounds, we have 
1 1 



■P(^t > x) -EYov[x,oo) 



^ v\x^ oo) 




^ETt - EYo 
t 


t 



< 



v\x^ oo) 



^ \EY,-EYo\ds+^E{Ti 



Next, (|26p and Jensen's inequality imply that 

1 /■* 1 1 /■* 

limsup — / \EYs - EFqI ds < oo, limsup -irE{Tl) < limsup - / EY^ds < oo, 

t->-o t Jq t_j.o t t^o t Jo 
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and (|28|) will follow. In order to show (|29p . consider now 

1 f 1 



t t 



d2ix) 
2 ' 



and note that, in view of Theorem 12.11 there exist constants to{e) > and K E (0,oo) such that 



sup|G,(t)| <i^, and |G,(t)|<e, 



for any < t < to. As before, 



1 



1 



1 



f2 



i2 ^ i2 



1 



Tt 
d2{x) 
2t2 



z/[x, CX)) 
t2 



E(T2) 



oo) 



The first term in the last expression can be bounded as follows: 
1 



< 



pd2{x) ^iy[x,oo) 



Then, it is now clear that we can bound the expression 



^V[Zt>x)-jEYov[x,^) 



as follows 



A < e^E{T^) + K^E{T^l{T,>t,}) + y{x, oo) 



'^2 y ti ' 12 

The third term on the right hand side of the above inequality is such that 



\d2{x)\ 



^E(T2) - p 



1 /I 



t \ t 



-ETt - EFo 



2 = i 
2 f2 



si -{EYs - EYo) -7 Ws, 



which converges to as t —t' due to (iii) in ([26|) . Hence, using (iv)-(v) in ([27]) and 



lEm'l{T,>to}) < E(T3)/to < e / E(r/)ds/io, 



we have 



limsup-Dt < ep, 
t->-o 



which implies ([^^ because e is arbitrary. ■ 

Proof of Lemma l4.3l We start by introducing some notation. Suppose that, under P* , X has Levy-Ito decomposition 



Xt = b*t 



z p*{dz, ds) 



J\z\<l 



z p* {dz, ds), 



(A-5) 



J\z\>l 



where p* is an independent Poisson measure on R\{0} x M-|- with mean measure i'* (dz)dt, and p* (dz, dt) := p* {dz, dt) - 
v*{dz)dt. Next, for a given fixed £ > 0, we set 



'^{\z\>e}tJ'*{dz,ds), and X^ := Xt - X^ 



(A-6) 



23 



hence, X"^ is a compound Poisson process with intensity :— v*{\z\ > e) and jumps {^f}; with common distribution 
l|2|>£i'*((iz)/Ae, while the remainder process X"^ is a Levy process with triplet (0, 6*, l{|2|<e}j^*((iz)), where 



b; := b* 

Let us fixed e = x/2. We first note that 



/ z'i-{\z\>6}''*{dz). 

J\z\<l 



for any t, a; > and for some universal constant K. Indeed, if we let denote the number of jumps before time t of 
the compound Poisson process X^, then we have 

P* [Xl > < P* {Nf 7^ 0) = 1 - e-^'* < K_t = v{{z : \z\ > x/2})t < Cx'^t. 

We now estimate P* (Xf > x). First, note that, due to the symmetry of the Levy measure v, 



E*iX!)^tibl+ zl{|,|<,}Z.*(dz)) = t(fe*- / zl{|,|>,}e^:/(dz)+ / zl{|,|<,}e^i/(dz)) 

■ ■ "'NI<1 •^|2|>l 





- / ^i{i 




J\z\>l 


t{b+ 






J\z\<l 


t{b+ 






J\z\<e 



t{b+ / z{e'' - l)v{dz) - I zl{|2|>5}e''i'(dz) + / zli^\^\<^}e''v{dz)) 
■ ■ 'kl<i •^|z|>l 

z(e" - l)v{dz)). 

\z\<x/2 

Thus, using concentration inequalities for centered random variable (e.g. }Hou02j . Corollary 1), for x > 2EXf , 



9p ( 9k ~^ 



2 



(Xf > a;) < P*(X; - E*Xf > x/2) < e"^ ' ; "°V* ' ='^^7 < ( ni^ ) < ^^i, 

\ ex J X'' 

where := Var*(Xf) = J^^^^^^y z'^v*{dz). Since M > 1, there exists a universal constant K such that 

2 9 2 

2C/;/^^^^^dz _ 2C(x/2)^-^ _ 

x^ i2-Y)x^ ■ 

We conclude that P*(Xf > x) < Ktx~^ for i(6 + /|^|<j./2 ^(^"^ — ^)'^{dz)) < x/2. This completes the proof, since 

F*{Xt >x)< P*{X^ > x/2)+¥*{X^ >x/2)< Kix'^ , 
whenever t{b + /|^|<j,/4 ^(e^ — l)i'{dz)) < x/4. ■ 

Proof of Proposition l472l Without loss of generality, we assume 5*0 = 1. We break the proof into two parts: 

(1) Let us assume through this part that {Xt)t is a symmetric CGMY process. Let b{u) := 5 + /|2|<„-2^(6^ ~ l)i^{dz). 
Obviously, 



5(u) < |6(it)| < |6| + / \z\\e' - l\iy{dz) + 2 \z\e' vidz) -.^^ b < oo. (A-7) 

J\z\<i A^\>^ 

Next, we write 

oo poo 

e"*''''"P* (Xt > ii/^uj = J l{„/4<ti^v^fa}e"*'^''"P* (X* > t'/^^w) (A-8) 

+ / l{«/4>*i-i/-5}e"*''""P* (Xt > t'/'^u) du. (A-9) 
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Clearly, e *^''''"P* (^Xt > t^^^u) < 1, so the first term converges to as t — ;> because Y e (1,2). From the inequality 
(|A-7p . we have 

l{ti/4>ti-i/>'6} = ^{t^/Yu>tb} - -'-{ti/>'„/4>t6(ti/y„/4)}' 

and using Lemma 14.31 we obtain that 

< mm{K{t^^^uy^t, 1} = min{A'w^'*', 1}, 

which is intcgrable because Y e (1,2). Hence, we can apply dominated convergence in the second term (|A-9|) and, 
using Proposition 14. 11 wc obtain that 



lim 



'{Z > u)du = E*(Z+). 



This show the result in view of ([Mt - pO)) . 

(2) In this second part, we relax the symmetry restriction. The idea is to reduce the problem to the symmetric case 
by applying a change of probability measure0 Concretely, let /3 := and, as in the proof of Proposition [521 define 

a probability measure P on (fi, J^) such that 



\b) =E(e'3-^''lB) /E (e''^*) , 



(A-10) 



for any B £ Ft- We can check that, under P, {Xt)t is a symmetric CGMY model with C — C ^ Y ^ Y , and 
G^M = {M + G)/2. Indeed, it follows that 

E (e™^*) = E (ef^+W^'j /E (e-^^') 

= exp \t CT(-y) {{M - 13-iuY + {G + 13 + iu)^ - (A/ - /3)^ - (G + /3)^' } + ibut 

Also, assuming /? > 0, 



E 



((e^' - 1)^ - E (e^' - l) + | = E ((e^' - l)^ (e^^' - l)) < ^ ((e^' - l) (e^^' - l)) = 0(t), 



since the moment function 1^9(2;) :— j{e^ — l)(e^^ — 1) ^ x"^ and Theorem l.l-(ii) in |FL08j can be applied. If /3 < 0, 
then 



E 



((e^' - 1)^ e'^^*) - E (e^'' - l) + | = E ((e^' - l)+ (l - e'^^*)) < E ((e^* - l) (l - e^^'')) = 0(i), 



for the same reason. Then, we only need to consider the asymptotic behavior of E (^(e^* — l)^ e^^*j as t 0, because 
Y € (1,2) so the 0{t) terms above are smaller than 0{t^^^). However, 



E 



1)- 



E I 



E( (e^* -1) 



and thus, using the fact that {Xt)t is symmetric under P and part (1) in this proof. 



limi'i/>'E((e^' - l)+e'^^"') = J™ ^^'^''^ ((e""* " ^) +) =^*(^- 



similar argument is applied in the proof of Proposition 5-(2) in |TnkvlO] but with a different aim. 
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